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The reversible Frank model is capable of amplifying the initial small statistical deviations from 
the idealized racemic composition. This temporary amplification can be interpreted as a chiral 
excursion in a dynamic phase space. It is well known that if the system is open to matter and 
energy exchange, a permanently chiral state can be reached asymptotically, while the final state is 
necessarily racemic if the system is closed. In this work, we combine phase space analysis, stability 
analysis and numerical simulations to study the initial chiral excursions and determine how they 
depend on whether the system is open, semi-open or closed. 
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I. INTRODUCTION 

The Prank model [l[ has been extensively invoked to justify theoretically the emergence of biological homochiral- 
ity @, H[ , and is usually analyzed as a reaction network in open systems (matter and energy are exchanged with the 
surroundings) composed of an irreversible enantioselective autocatalysis coupled to an irreversible mutual inhibition 
reaction between the product enantiomers. The model shows how homochirality is achieved as a stationary state 
when the mutual inhibition product (the heterodimer) is removed from the system and when the concentration of the 
achiral substrate is held constant. By contrast, for reversible transformations and when the mutual inhibition product 
remains in the system, the final stable state can only be the racemic one. As a consequence, a thermodynamically 
controlled spontaneous mirror symmetry breaking (SMSB) cannot be expected to take place. In particular, SMSB is 
not expected for reversible reactions taking place in systems closed to matter and energy flow. 

Nevertheless, as was recently demonstrated [|J for systems closed to matter flow, the Frank model is a prime 
candidate for the fundamental reaction network necessary for reproducing the key experimental features reported on 
absolute asymmetric synthesis in the absence of any chiral polarization [5] . Most importantly, when reversible steps 
in all the reactions are allowed it is capable of [ij (i) amplification of the initially tiny statistical enantiomeric excesses 
from ee ~ 10 _8 % to practically 100%, leading to (ii) long duration chiral excursions or chiral pulses away from the 
racemic state at nearly 100% ee, followed by, (iii) the final approach to the stable racemic state for which ee = 0, 
i.e., mirror symmetry is recovered permanently. To understand this temporary asymmetric amplification is important 
because the racemization time scale can be much longer than that for the complete conversion of the achiral substrate 
into enantiomers. 

Long duration chiral excursions have also been reported recently in closed chiral polymerization models with re- 
versible reactions [(| where constraints implied by microreversibility have been taken into account. These results 
are important because they suggest that temporary spontaneous mirror symmetry breaking in experimental chiral 
polymerization can take place, and with observable and large chiral excesses without the need to introduce chiral 
initiators Q or large initial chiral excesses Q. 

The purpose of this Letter is to elucidate the nature of these chiral excursions by combining the information provided 
by phase plane portraits, numerical simulation and linear stability analysis. We consider the Frank model, this being 
the most amenable to such types of analysis and because it is the "common denominator" of numerous more elaborate 
theoretical models of SMSB Q. 

The reaction scheme consists of a straight non-catalyzed reaction Eq. ([TJ , an enantioselective autocatalysis Eq. © , 
where A is a prechiral starting product, and L and D are the two enantiomers of the chiral product. We also assume 
reversible heterodimerization step in Eq. where LD is the achiral heterodimer. The ki denote the reaction rate 
constants. In the following, we give the reaction steps in detail. 
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Production of chiral compound: 



Autocatalytic amplification: 



k\ k\ 



A Z L > A Z D - (1) 



k 2 k 2 

L + A,^L + L, D + A^D + D. (2) 

Hetero-dimerization: 

L + D "~ LD. (3) 

k- 5 

We assume the feasibility of the reverse reaction for all the steps. Focusing our attention on chiral excursions, we 
make a careful distinction between open, semi-open or fully closed systems. These system constraints are crucial for 
determining both the intermediate and the asymptotic final states of the chemical system. 



II. OPEN SYSTEM 
A. Rate equations 

We first consider the original Frank scenario There, steady and stable chiral states can be achieved, since the 
system is permanently held out of equilibrium. See [i[ for more details. An important question is, can the system 
support chiral excursions? That is, pass through temporary chiral states before ending up in the final racemic state? 
In the original Frank model there is an incoming flow of achiral compound A and elimination of the heterodimer LD 
from the system. A convenient way to account for the inflow of achiral matter is to assume that the concentration 
of the prechiral component [A] is constant, and then we need not write the corresponding kinetic equation for it. 
For the outflow the heterodimer leaves the system at a rate j. We assume that the heterodimer formation step is 
irreversible, and set fc_s = 0. Note, the elimination of LD from the system can actually be neglected as long as the 
hetero-dimerization step is irreversible Q. We retain this outflow however since it is needed to obtain stationary 
asymptotic values of all three concentrations [L], [D] and [LD]; see the fixed points below. So with [A] = const and 
replacing Eq. §5\ ) by 

L + D H LD, (4) 
LD 4 0, (5) 



we obtain the rate equations 



j t [L] = h[A] + {k 2 [A] - k-!)[L] - k_ 2 [L] 2 - k 5 [L][D], (6) 

j\\D] = k t [A] + (k 2 [A] - k-i)[D] - k_ 2 [D] 2 - k 5 [D][L], (7) 

j t [LD] = k 5 [L][D]-j[LD]. (8) 
The key variable throughout is the chiral polarization 

V [L] + [D] ' (9j 

also called enantiomeric excess ee, which obeys —1 < n < 1 and which represents the order parameter for mirror 
symmetry breaking. 

In order to simplify the analysis, we define a dimensionless time parameter t = (fc2[^4] — and dimensionless 

concentrations that scale as [L] = fc 5 (fc 2 [^4] - fc-i) -1 ^ [D] = k 5 (k 2 [A] - fc_i) _1 [L>], [LD] = k 5 (k 2 [A] - fc_i) _1 [L£)]. 
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It is convenient to define the sums and differences of concentrations: \ — [L] + [D], y = [L] — [D], and for the 
heterodimer put P = [LD\. The chiral polarization r\ — yjx remains unchanged. 
In terms of the new variables, Eqs. ([6][8]) read 

% = 2u + X-i(.9 + l)x 2 -^(.9-l)xV, (10) 
dP 1 o, 

to = - 4 X\1-V 2 )-7P- (12) 
The dimensionless parameters appearing here are: 

kik 5 [A] fc_ 2 7 

f> 9 = -T- > 7 = 777771 il \ • ( 13 ) 



(fc 2 [A]-fc_i)2' a fc 5 ' ' (fe 2 L4] - fc_0 

The system is described by three equations Eqs. (|10H12p . Since P does not enter into the equations for \ ar± d f], the 
equations decouple and the dynamical system to study is effectively two-dimensional and so the appearance of SMSB 
cannot depend on whether the heterodimer is removed from the system when fc_5 = 0, although the fixed points will 
(see Sec ETC]). 

B. Phase plane and linear stability analysis 

In the phase space of the dynamical system defined by Eqs. ([lOllip there are curves with a special significance. 
These are the nullclines defined by 

fr = °- (15) 

The intersections of these curves give the possible steady states (or fixed points) of the system. The condition 
dx/dr = leads to two curves 

(1) l±yT+ Mg(l + r?) + rrg) 
X± g(l + r) 2 ) + l-r] 2 ' 1 J 

while drj/dr — implies the three curves 



V 1 - 9 + (9 - W 

For it > the solutions denoted X- correspond to negative total enantiomer concentrations so we discard them. 
The physically acceptable nullclines are plotted in Fig. [T] Which of the two different intersection configurations is 
obtained depends only on the single parameter g. We emphasize that despite the similar appearance, the nullcline 
graphs should not be confused with the classic bifurcation diagrams that have been discussed often in the past [Tol - [T^ | . 

C. Fixed Points and Stability 

The system has several steady states: besides an unphysical state that we disregard, there is a Z 2 pair of chiral 
solutions Q±, and a racemic state R: 



4it 



(p 2( g + l)u+ A /4( g +l)n + l + l 1 + y/4(g + l)u+l \ 

R = [ p = oi _l >x= — - — — j ,y = o), (is) 

2(5 + 1)^7 g+1 I 

Q± = [P= ,x = -,y = ± i=^r ■ (19 

7-57 9 V9 - 1 ' 
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FIG. 1: Nullclines for the open system (| 101 11 p . These curves correspond to u = 0.3 implying g cr it = 0.59. The dr\jdr = 
and dx/dr = nullclines are plotted in black and red, respectively. The black (red) arrows indicate the regions of phase space 
where r\ (x) increases or decreases. Left: g = 0.79 > g cr it- The nullclines intersect in only the one point R representing the 
asymptotic stable racemic state. Right: g — 0.29 < g cr it- In this case there are three intersections, Q± and R representing the 
Z% equivalent stable chiral states and the unstable racemic state, respectively. From Eq. (|22[) the enantiomeric excess at Q± is 
77 = ±0.93. 



The associated eigenvalues are given by @: 



W*) = (-,/ 4( fl+ l) U + l , 1 - g ^ 1)u+1 ,-7), (20) 

^ -y/l6ff 3 M + V -4.9 + 1-1 ^16ff% + 4g 2 - 4 g + 1 - 1 \ 
Ai, 3 , 8 (Q±) = ( Yg , Yg -l). (21) 

Note that \\{R) < and Ai(Q) < are always negative whereas A2 (ii) > and Aa(Q) < for 5 < 5 cr i t , otherwise 
A2(i?) < and A2(Q) > for g > g cr a, where g C rit = (VI + 16u — l)/8u is the critical value for this parameter. Note 
that g C rit{u) < 1 for all u > 0. For small u we can write g cr it — 1 — 4m; while for large u, g C rit —> 1/2U 1 / 2 . Thus the 
direct chiral monomer production step (oc k± in ([!])) tends to racemize the system leading to final 77 values strictly 
less than unity: 



which holds when g < g cr n. The chiral monomer production step thus reduces the range of g for which stable mirror 
symmetry breaking can occur, and the chiral solutions are no longer 100% chiral. 

Figure [5] shows the temporal evolution of the L and D chiral monomers starting from an extremely dilute total 
enantiomer concentration and the very small statistical chiral deviations from the ideal racemic composition. The 
right hand side of this figure shows the evolution in terms of the quantities x an d V- Note the mirror symmetry 
breaking signalled by r\. Are there chiral excursions found in the open system model? A chiral excursion holds when 
the enantiomeric excess 77 departs from a small initial value, evolves to some maximum absolute value and then decays 
to the final value of zero. To ensure a final racemic state we must set g > g cr it but then we find no numerical evidence 
for such temporary chiral excursions. This can be understood qualitatively from inspection of left hand side of Fig. [T] 
The initial conditions (dilute chiral monomer concentration and statistical chiral fluctuation) corresponds to a initial 
point located at tiny values of x an d close to the vertical nullcline, well below the point labeled as R. The system 
is attracted to the black curve and moves up the curve to R. In this situation, it is impossible for the chiral excess 
to increase, not even temporarily. Notice the time scales for x an d 77 are of the same order. On the other hand, if 
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FIG. 2: Chiral symmetry breaking in the open system ([6][8]). Temporal behavior (logarithmic scale) of the individual enantiomer 
concentrations [L] and [D] (left), and the chiral polarization r\ and total enantiomer concentration \ (right)- Initial concentra- 
tions: [L] = (1 x 10~ 6 -fix 10" 15 )A/, [D] = 1 x 10" 6 M (770 = 5 x 10~ 8 %) and [A] = 1M. Rate constants: fei = 10" 4 s" " 



k-i = lO^s" 1 , k 2 = ls^M' 1 , k-2 = O^s^M -1 and k s = ls _1 Af _1 . 
figures of simulations, we always display original variables [L], [D], [LD] 



These values correspond to g = 0.5 and u — 10 
etc. as function of time t. 



In 



9 < 9crit, then we have the situation depicted on the right hand side of the figure. Here the same initial point moves 
towards the vertical nullcline and up towards R, but once past the locally horizontal black curve, is attracted to one 
of the two chiral fixed points where it stays forever, provided the system is maintained out of equilibrium. The chiral 
symmetry is permanently broken, and there is no excursion such as we have defined it. 

III. SEMI-OPEN SYSTEM 
A. Rate equations 

To elucidate the temporal evolution of x an d V f° r a more general setting, we do not remove the heterodimer from 
the system now allow the back reaction to chiral monomers, and we keep [A] constant. There is an implicit inflow as 
a consequence of constant [A], but no outflow, so we denote this case "semi-open". While there is no mass balance 
the system can still exhibit temporary SMSB. 

The corresponding rate equations are 

^-[L] = h[A] + (k 2 [A] - k-{)[L] - k^ 2 [L} 2 - k 5 [L][D] + fc_ 5 [££>], (23) 
at 

j t [D] = k 1 [A] + {k 2 [A]-k. l )[D]-k- 2 [D} 2 - k 5 [D][L] + k_ 5 [LD], (24) 

1[LD] = h[L][D}-k- 5 [LD]. (25) 
at 

After performing the same rescaling as in the open case, we arrive at 

f: = 2^ + x - \(g + i)x 2 -\{g- i)xV + 2rP, (26) 

(27) 

' !P \ 2 (l-V 2 )-rP. (28) 



dr 4 

The dimensionless parameters appearing here are: 

hh[A] 



>9 = -?;r= fuI -\ r . (29) 



(k 2 [A] - k_,)^ " k 5 ' (k 2 [A\ - fe_i) 



The system is described by three equations Eqs. (I26H28P which do not decouple. 
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B. Phase plane and linear stability analysis 



In order to obtain an approximate two-dimensional phase plane representation of the system, we will invoke the 
dynamic steady state approximation for the heterodimer. Such approximations are usually justified when there 
exists a clear separation of time scales in the problem, thus allowing one to identify rapidly and slowly changing 
concentrations |l3j. Here however, no such time scales are evident, all concentration variables evolve on a similar 
time scale. Nevertheless, we will see a posteriori that this approximation can be good over a wide range of time 
scales. We therefore assume that the heterodimer is in a dynamic steady state P s tat relative to the chiral monomer 
concentrations and chiral polarization: 

Pstat* ^i 1 -^)- (30) 

Substituting this P s t a t into Eqs. (I26I27[) leads to the differential equations 

^ = 2u + X - 9 - X i(l + r,% (31) 



air ~ 2 

| =<-«-7 + fxd + ^)). (32) 



As above, we study the phase space of the two-dimensional system by means of the nullclines. These are plotted in 
Fig. [3l The condition dx/dr — implies two curves 



l±y/T+ 4gu(l + ^) 

X± = TT, — T\ ' ( 33 

9(1 + T) 



whereas the condition drj/dr — implies the three curves 



7 ? = 0, ^ ± =±/l+i!i. (34) 

9X 




FIG. 3: Nullclines for the semi-open case (|31I32[) in the steady state approximation for P. The nullclines (|33I34[1 are plotted 
in red and black, respectively. Red and black arrows indicate the phase-space regions of increasing or decreasing \ an d V- The 
four stationary solutions O, R, Q± are indicated by the black dots. These curves are illustrated for u = and g — 0.5. 
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C. Fixed Points and Stability 



We solve Eqs. (|26fl28[) looking for steady states. To keep the algebra manageable, we also set u = as in [9(. There 
are four solutions, namely, the empty O solution, the racemic R or the two mirror-symmetric chiral Q± solutions: 

O = (P = 0,x = 0,I/ = 0), (35) 

R = (p=-^-,X=-,V = o) , (36) 
V 9 r g J 

Q± = (p = 0, x = -,y = ±-). (37) 

Note that the final heterodimer concentration P is zero in the chiral states Q±. Note that the steady state approx- 
imation Eq. ([50)1 implies the same result since \r]\ = 1. In order to study the linear stability of the four possible 
homogeneous solutions 0,R and Q±, we calculate the eigenvalues of the 3x3 Jacobian array M open derived in @ 
after deleting the 3rd and 4th rows and columns. The eigenvalues corresponding to these solutions are given by 

Ai, 2 , 3 (0) = (1,1, -r), (38) 



/ 2 + g(l + r) + y/4 + g 2 (-l + r) 2 + 4g(l + r) 
Ai Aa (J2) = (-1, Yq ' 



-2 - g(l + r) + y/4 + ff 2 (-l + r) 2 + 4g(l + r) 
2^ 



2g 

"l 2 4- drill +r\\ 

(39) 



Ai,2,a(Q±) = 



/_ 1 + g (-l + r) + y/l + 2 g (-l + r) + g 2 (l + r) 
^ ~ V ' 2o 



2g 

-1 + ,g(l - r) + V(l + g(-l + r)) 2 + 4g 2 r 
25 



). (40) 



As Ai_2(0) > 0, the empty state is always unstable. An inequality analysis shows that both \2{R) < and \s{R) < 
for all r > and g > 0. Since Xi(R) = — 1 this demonstrates that the racemic state R is always stable. As an 
independent check, we also verify that As(Q) > is positive for all r > and g > 0, so the chiral solutions Q± are 
always unstable. The final outcome will always the racemic state. There is no stable mirror symmetry broken solution 
when the heterodimer dissociates back into the chiral monomers. Nevertheless, the system can have temporary chiral 
excursions. 




FIG. 4: Temporary chiral symmetry breaking and chiral excursions in the semi-open system (|23I25|) . Temporal behavior 
(logarithmic scale) of the individual enantiomer concentrations [L] and [D] (left) and the chiral polarization r\ and total 
enantiomer concentration \ (right). Initial concentrations: [L]o = (1 x 1CF 7 + 1 x 1CF 15 )M, [D]q = 1x W~ 7 M (r/o = 5 x 10 _8 %) 
and [A] = 1M. Rate constants: fei = 10" 4 s"\ k-i = 10" 6 s _1 , fc 2 = ls^M -1 , fc_ 2 = 0.5s _1 M _1 , k 5 = ls~ l M~ l and 
ks — 10 _3 s _1 . These rate constants imply g = 0.5 and u = 10~ 4 . 



In Fig. @] we plot the temporal evolution of the L and D chiral monomers starting from an extremely dilute total 
enantiomer concentration and the very small statistical chiral deviations from the ideal racemic composition. The 
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FIG. 5: Semi-open system: Comparison of direct numerical solution P = [LD] and the steady state approximation P s tat in 
Eq. (|30[) for the heterodimer concentration (after transforming from [L], [D], [LD] to Xi^-P)- 



right hand side of this figure shows the evolution in terms of the quantities \ and rj. Note the chiral excursion in rj 
for the time interval between t ~ 100s and t ~ 1000s. Finally, we compare the heterodimer concentration from direct 
numerical simulation with the steady state approximation in Fig. [5] To do so we simulate the original concentration 
variables in Eqs. (|23H25l) and then transform results in terms of %, rj and P. While P s tat appears to overestimate 
P = [LD], it provides a reasonably good approximation to the actual heterodimer concentration P right after chiral 
symmetry is broken, at about t ~ 100s and coincides perfectly after chiral symmetry is recovered, and when x reaches 
its asymptotic value, after approximately t ~ 10 4 s. 



IV. CLOSED SYSTEM 
A. Rate equations 

The rate equations directly follow from (p][3]): 

j t [L] = h[A] + [k 2 [A] - k-i)[L] - k^ 2 [L] 2 - k 5 [L][D] + k- B [LD], (41) 

-J [D] = h [A] + (k 2 [A]-k_ 1 ) [D] -k_ 2 [D} 2 - k 5 [D] [L] + fc_ 5 [LD] , (42) 
at 

j t [A] = -2k x [A] (k 2 [A] - fc_0([i] + [D]) + fc_ 2 ([L] 2 + [Df), (43) 

j\\LD] - k 5 [L][D]-k- 5 [LD]. (44) 

There is no flow of material into or out of the system. Since [A] is not constant in this situation, we cannot use it to 
rescale the time or the concentrations. Instead, we take r = kit for the time parameter and [L] = (fcg/fci)[i], etc. for 
the dimensionless concentrations. This allows us to express the rate equations in the following dimensionless form: 

±[L] = [A]-u[L]+h[A][L]-g[L] 2 -[L][b]+p[LD], (45) 
4-[D] = [A] - u[D] + h[A][D] - g[D] 2 - [D][L] + p[LD], (46) 

(IT ... 

±[LD] = [L][D]-p[LD]. (47) 

These are subject to the constraint [A] = [C] — ([L] + [D]) — 2[LD], where [C] = (kc,/ki)[C] and [C] is the total initial 
concentration, being constant in time. The four parameters appearing here are 

k-i k-2 , k 2 k- 5 

t\i± n<5 ft/5 ">1 
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Changing variables as before to x, Vi Pi we arrive at 



^ = 2[I] + (h[A] - u) x - \{g + i)x 2 -\{g- i)xV + 2pP, (49) 
= v (h[A]-g X -u)-l(Q), (50) 



dr ' x \dr 

f = \x 2 (l-V 2 )- P P (51) 
In these variables, the constant mass constraint reads [A] = [C] — x — 2P. 

B. Phase plane and linear stability analysis 

As in the semi-open case, to obtain an approximate two-dimensional phase plane portrait, we assume that the 
heterodimer is in an approximate steady state P s tat and solve Eq. (|51J for 

P sta t~^(l-V 2 )- (52) 
Substituting this back into Eqs. (|49l50p . we obtain 

^ = (2 + h X )[A] -ux- |x 2 (l + v 2 )i (53) 
dv ( 2[1] g 

where 



X(l + ^ 2 )), (54) 



[A] = [q-x-^(i-^ 2 )- (55) 

The nullcline condition dx/dr — leads to an unwieldy cubic equation in x- More importantly, the nullcline is an 
even function of 77, reflecting the underlying Z2 mirror symmetry. The other condition drj/dr = is straightforward 
to solve analytically and leads - after discarding the unphysical solution corresponding to negative total enantiomer 
concentrations - to two curves 

_ -1 + + 2[C]( g (l + jgj/2 + (1 - f)/P - 9) 
V ' X_ g(l + r?)/2 + (l-r,*)/p-g ■ { > 

The curve xiv) is an even function of 77 and xi 7 ]) ~~ ^ 00 f° r V 2 ~ * 1 with a minimum at 7/ = 0. Thus the nullcline has 
the form of a narrow fork as depicted in Fig. [6] 



C. Fixed Points and Stability 



A linear stability analysis for the closed Frank model is given in [9J. That analysis was carried out in terms of Xi y 
and P and does not assume the stationary approximation for the heterodimer. It turns out, even in a model as simple 
as this one, that keeping p > is analytically untractable, so we consider p = in what follows. Then the asymptotic 
stationary racemic R and chiral solutions Q± are given by 

R = (p = ^Y-,x = o,y = o), (57) 

„ („ \C]q + u u u\ , , 

Q± = [P= LJ % ,X=-~,V = ±-)- 58 

v 2 ff g gJ 

The chiral solutions Q± are unphysical for all u > since they imply negative total enantiomer concentrations x < 0. 
Thus, the only physically acceptable solution is the racemic one R, and this is (at least marginally) stable, the 
corresponding eigenvalue was calculated to be [9j 



Ai, 2 , 3 (i?) = (0,-2 -«,-«). 



(59) 
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FIG. 6: Nullclines for the closed system (|53l54p in the steady state approximation for P. The r\ and \ nullclines are plotted as 
the black and red curve, respectively. These nullclines intersect at the one point (indicated with a dot) which corresponds to 
the stable racemic solution R with r\ — and x > 0. The flow directions are indicated with arrows. This diagram corresponds 
to the parameter values g = 0.05, u — 0.01, h — 0.1 and p = 1.0 and [C] = 10°. 
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FIG. 7: Temporary chiral symmetry breaking and chiral excursions in the closed system (|41M2p . Temporal behavior (loga- 
rithmic scale) of the individual enantiomer concentrations [L] and [D\ (left) and the chiral polarization r\ and total enantiomer 
concentration \ (right). Initial concentrations: [L] = (1 X 10" 7 + 1 X 10~ 15 )M, [D] = 1 x 10" 7 M (r] = 5 x 10" 8 %) 
and [A]o = 1M. Rate constants: ki = 10~ 4 s~\ k-i = 10 _6 s _1 , k 2 = ls _1 M _1 , fc_ 2 = 0.5s' 1 M" 1 , k 5 = lOs^M" 1 and 
fc_5 = 10- 4 s _1 . 



In the limit p = 0, the substrate is consumed and all the matter ends up finally as pure heterodimer. Finally, note that 
lim u _>.o Q± = R- the unphysical chiral solutions merge to the racemic one when = 0. For reversible heterodimer, 
the matter in the racemic state 77 = is distributed between the chiral monomers and the heterodimer: P = ([C] — x) /2 
and x > in keeping with the law of mass action. The single intersection R displayed in Fig. [6] indicates that the 
racemic state is the only possible solution, in qualitative agreement with the stability analysis. 

Example of a chiral excursion in a closed system is provided in Fig. [7] Once again, as for the semi-open situation, 
the scheme is capable of amplifying a tiny initial chiral excess to practically 100%, followed by final approach to the 
racemic state. The steady state approximation for the heterodimer is rather poor during the early stages (Fig. [8]), but 
similar to semi-open case, converges to the true heterodimer concentration after symmetry breaking and restoration. 
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FIG. 8: Comparison of direct numerical solution [LD] and the steady state approximation [LD] ata t = (fc/fc-s )[£■][!?] for the 
heterodimer concentration. The steady state approximation peaks to a maximum value of about 25 in this example. 



In this Letter, we investigated transient mirror symmetry breaking in chiral systems, in particular in the Frank 
model in settings of open, semi-open, and closed environments. Temporary chiral excursions are observed for closed 
and semi-open systems and explained through phase space analysis, stability analysis and numerical simulations. 
Such chiral excursions may be experimentally observed and could be mistaken for a transition to a chiral state. They 
are in fact a long sought goal for the experimental chemist who could actually fail to see them if not aware of their 
transitory nature. In open systems by contrast, the racemic state is approached monotonically. Therefore, it is 
important to understand the processes and constraints responsible for these outcomes. This Letter has focused on 
the effects that the in- and outflow of matter has on these phenomena. The open nature of the Frank model can 
be arranged experimentally with an incoming flow of achiral precursor A and an outflow of the product heterodimer 
LD, to conserve mass balance. Mathematically, we can model the inflow by assuming a constant concentration of [A] 
and the outflow by a term representing the rate at which LD leaves the system. In our open model we assumed no 
dissociation of LD back into chiral monomers. From the point of view of achieving permanent SMSB, there is then 
actually no need to remove LD from the system, but we retain this outflow since it is needed to ensure stationary fixed 
points for all the chemical concentrations. For the semi-open case, LD is not removed and we allow for its dissociation 
into chiral monomers. There is no mass balance but temporary symmetry breaking can arise. Finally, in the closed 
system there is neither inflow nor outflow, total mass is conserved and temporary symmetry breaking can occur. 

A recent kinetic analysis of the Frank model in closed systems applied to the Soai reaction Q indicates that in 
actual chemical scenarios, reaction networks that exhibit SMSB are extremely sensitive to chiral inductions due to 
the presence of inherent tiny initial enantiomeric excesses (J|. This amplification feature is also operative in much 
more involved reaction networks such as chiral polymerization When the system is subject to a very small 
perturbation about an extremely dilute racemic state, the initial chiral fluctuation does not immediately decay, but 
becomes amplified and drives the system along a long-lived chiral excursion in phase space before final and inevitable 
approach to the stable racemic solution. Mauksch and Tsogoeva have also previously indicated that chirality could 
appear as the result of a temporary asymmetric amplification fbl . [l5j . 

Excursions in phase space as studied in this work are superficially reminiscent of excitable systems as studied in 
dynamical systems [HI 0, [ljj ■ But there are important differences. First of all, the excursions reported in chiral 
systems are not easily visualized in the chiral monomer concentrations themselves, but are strikingly manifested by 
the chiral polarization or enantiomeric excess. Secondly, the total enantiomer and heterodimer concentrations do 
increase with time, so that the complete phase-space trajectory does not follow a closed path: there is no return to 
the initial state. The chiral excursion is a one way trip, not a round trip as in an excitable system. Third, whereas 
excursions are traditionally studied for open excitable systems [TH, [TH, [I?} , chiral excursions are observed here only 
for closed or at most semi-open systems, but not for open systems. 

The original impetus for considering phase-space descriptions of the Frank model comes not only from the chiral 
excursions reported in Q and [|| but also by the recent report of damped chiral oscillations detected numerically 
in a model of chiral polymerization in closed systems [f|. Absolute asymmetric synthesis is achieved in the latter 
scheme, accompanied by long duration chiral excursions in the enantiomeric excesses for all the homopolymer chains 
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formed, analogous to the much simpler Frank model. But unlike the latter, strong enantiomeric inhibition converts 
these excursions into long period damped chiral oscillations in the enantiomeric excesses associated with the longest 
homochiral polymer chains formed. Moreover, short period sustained chiral oscillations have been observed numerically 
in a recycled Frank model open to energy flow, for large values of the inhibition [l8[ . This oscillatory behavior poses 
an additional problem for the origin of biological homochirality, since any memory of the sign of the initial fluctuation 
is further erased by subsequent oscillations thus adding a further element of uncertainty to the overall problem. 
Chemical oscillations have been traditionally studied in conjunction with excitability. Although the latter concept is 
not directly applicable to models exhibiting SMSB, it remains to be seen if the techniques used to study oscillations 
can be applied profitably to reaction schemes that lead to chiral oscillations. 
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